REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of 

gethering  and  melnteining  the  date  needed,  jnd  c< 
of  informetion,  including  suggestions  for  reduc 
(0704-0188)  1215  Jefferson  Davis  Highway,  Su 
subject  to  eny  penalty  for  felling  to  comply  with  e 

PLEASE  DO  NOT  RETURN  YOUR  FOF 

information  is  estimated  to  average  1  hour  per  response.  Including  the  time  for  reviewing  instructions,  searching  existing  deta  sources 

ompleting  and  reviewing  the  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  eny  other  espect  of  this  collection 
ing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Informetion  Operations  and  Reports 
ite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  thet  notwithstanding  eny  other  provision  of  iew  no  person  shaii  be 
collection  of  information  if  It  does  not  dispiey  a  currently  velld  0MB  control  number 

m  TO  THE  A80VE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM- Y YYY) 

01012008 

2.  REPORT  TYPE 

Journal  Article 

3.  DATES  COVERED  (From  -  To} 

3D  underwater  imaging  using  vector  acoustic  sensors 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 

6I153N 


6.  AUTHOR(S) 

Dennis  Lindwall 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S}  AND  ADDRESS(ES) 

Naval  Research  Laboratory 
Marine  Geoacoustics  Division 
Stennis  Space  Center,  MS  39529 


8  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NRL/JA/743  0-06-2 


9.  SPONSORING/MONITORING  AGENCY  NAME{S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington  VA  22217-5000 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

ONR 


11  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 


5  0Oct°5',’|c" 


13.  SUPPLEMENTARY  NOTES 

GEOPHYSICS,  VOL  73,  NO.  1  (JANUARY-FEBRUARY  2008);  P-Ql-17,  11FIGS.  10.1190/1.2813347 


14.  ABSTRACT 


Marine  surveys  that  use  vector  acoustic  sensors  may  allow 
for  3D  imaging  of  underwater  environments  with  a  much 
smaller  amount  of  data  than  current  3D  hydrophone  surveys. 
Newly  developed  sensors  make  vector-acoustic-based  sur¬ 
veys  practical.  This  concept  is  demonstrated  with  data  from  a 
three-axis  accelerometer  and  a  collocated  hydrophone  in  an 
acoustic  water  tank  using  a  short-pulse  source  and  passive 
scattering  targets.  One  algorithm  rectifies  the  vector  data  with 
scalar  pressure  data  and  another  maps  the  vector  data  into  a 
3D  volume,  showing  several  slices  of  the  volume  images. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

A  dCTD apt 

IB.  NUMBER 

OF 

19a  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c  THIS  PAGE 

Ado  1 MAv 1 

Ur 

PAGES 

7 

Dennis  Lindwall 

Unclassified 

Unclassified 

Unclassified 

UU 

19b.  TELEPHONE  NUMBER  (Include  area  code} 

228-688-5306 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSi  Std.  Z39  18 


GEOPHYS1CS.VOL.73.NO.  1  (JANUARY-FEBRUARY  2008);  P.Q1  -Q7. 11  FIGS. 
10.1190/1.2813347 


3D  underwater  imaging  using  vector  acoustic  sensors 

Dennis  Lindwall1 


ABSTRACT 

Marine  surveys  that  use  vector  acoustic  sensors  may  allow 
for  3D  imaging  of  underwater  environments  with  a  much 
smaller  amount  of  data  than  current  3D  hydrophone  surveys. 
Newly  developed  sensors  make  vector-acoustic-based  sur¬ 
veys  practical.  This  concept  is  demonstrated  with  data  from  a 
three-axis  accelerometer  and  a  collocated  hydrophone  in  an 
acoustic  water  tank  using  a  short-pulse  source  and  passive 
scattering  targets.  One  algorithm  rectifies  the  vector  data  with 
scalar  pressure  data  and  another  maps  the  vector  data  into  a 
3D  volume,  showing  several  slices  of  the  volume  images. 
The  imaging  algorithm  maps  the  scattered  energy  using  the 
direction  and  traveltime  independently  for  each  source-re¬ 
ceiver  pair  rather  than  using  the  phase  coherence  methods 
common  in  exploration  seismology.  Imaging  a  more  com¬ 
plex  and  realistic  marine  environment  requires  vector  wave- 
field  decomposition  techniques  and  other  theoretical  devel¬ 
opments  but  may  allow  for  3D  vector-acoustic  seismic  sur¬ 
veys  using  logistics  similar  to  2D  surveys  that  use  conven¬ 
tional  hydrophones. 


INTRODUCTION 

Vector  acoustic  sensors  measure  either  the  particle  motion  or  the 
pressure  gradient  from  passing  acoustic  waves  and  produce  three-  or 
four-component  (3-C  or  4-C)  data  (Fahy,  1977,  1989;  Berliner  and 
Lindberg,  1996).  Vector  acoustic  data  include  directional  informa¬ 
tion  that  may  allow  for  3D  acoustic  imaging  with  significantly  less 
data  than  required  with  hydrophone  data.  Using  a  vector  sensor  with 
a  controlled-pulse  source,  one  can  determine  both  distance  and  di¬ 
rection  for  the  sound  source  as  well  as  scattering  and  reflection 
points. 

Geophones,  which  usually  measure  displacement,  have  been  used 
to  measure  local  particle  motion  in  the  solid  earth  since  the  start  of 
observational  seismology.  One  significant  technical  difficulty  with 
geophones  is  the  coupling  between  the  solid  earth  and  the  sensors  as 


well  as  the  device  resonances  (D’Spain  et  al.,  1991).  Using  geo¬ 
phones  in  water  poses  different  coupling  problems,  but  new  devel¬ 
opments  in  miniature  accelerometers  that  allow  for  the  construction 
of  small,  rigid,  neutrally  buoyant  sensors  allow  for  much  better  cou¬ 
pling  between  the  sensors  and  the  local  particle  motion  (Shipps  and 
Deng,  2003).  Towing  an  accelerometer-based  sensor  through  the 
water,  such  as  in  a  streamer,  would  require  a  more  practical  suspen¬ 
sion  solution  than  the  one  used  in  our  experiment  while  minimizing 
the  flow  noise.  The  three  separate  accelerometer  components,  as 
well  as  an  accompanying  hydrophone,  also  must  be  matched  in  their 
phase  and  frequency  response. 

Data  from  vector  acoustic  sensors  used  in  geoacoustic  surveys  of 
the  seafloor  provide  additional  directional  information  not  present  in 
scalar  pressure  data  from  traditional  hydrophones  (Lindwall,  2006). 
Hydrophone  arrays  for  use  in  seismic  profiling  and  seafloor  imaging 
surveys  such  as  side-scan  sonar  need  to  be  about  20  wavelengths 
long  to  provide  a  directional  resolution  of  1°  or  2°  (Kinsler  et  al., 
1982).  Unambiguous  determination  of  direction  requires  three  or¬ 
thogonal  arrays  of  scalar  sensor  arrays,  a  configuration  that  is  cum¬ 
bersome  and  impractical.  A  combined  vector  sensor  and  controlled 
source  give  both  distance  and  direction  for  the  sound  source  as  well 
as  scattering  and  reflection  points. 

In  principle,  a  single  sound  pulse  recorded  at  a  single  vector 
acoustic  sensor  can  provide  a  3D  acoustic  picture  of  the  local  envi¬ 
ronment.  In  practice,  numerous  ambiguities  are  caused  by  overlap¬ 
ping  arrivals  from  different  scattering  points.  These  ambiguities  can 
be  resolved  with  a  single  line  of  source  or  sensor  locations  in  most 
cases.  Data  from  a  linear  hydrophone  array  can,  at  best,  give  only  a 
2D  image  of  the  environment  from  the  same  survey  configuration. 
For  specular  reflections,  adequate  illumination  of  reflective  surfaces 
may  require  multiple  source  and  receiver  locations. 

In  this  paper,  1  show  how  a  vector  acoustic  sensor  may  be  used  for 
underwater  applications.  I  also  show  that  vector  acoustic  data  can  be 
acquired  with  existing  sensors  and  that  images  of  scatterers  in  a  3D 
volume  can  be  made  with  a  rather  small  amount  of  data.  This  is  veri¬ 
fied  by  an  experiment  using  an  accelerometer-based  vector  acoustic 
sensor  in  a  water  tank  with  a  short-pulse  source  and  passive  scatter¬ 
ing  targets.  I  present  an  algorithm  for  rectifying  the  vector  data  with 
scalar  pressure  data  and  an  algorithm  for  mapping  the  vector  data 
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into  a  3D  volume.  I  then  show  several  slices  of  the  resulting  volume 
images. 


THE  EXPERIMENT 

The  data  from  the  experiment  were  acquired  with  aTV-001  minia¬ 
ture  vector  sensor  designed  for  use  in  water  by  Roger  Richards  and 
others  of  the  NAVSEA  division  of  the  Naval  Underwater  Warfare 
Center  in  conjunction  with  Wilcoxon  Research  Inc.  (Shipps  and 
Deng,  2003).  The  sensor  contains  a  three-axis  accelerometer  and  a 
hydrophone,  has  a  sensitivity  of  1 .0  V /g  and  a  frequency  response  of 
about  3  Hz  to  9  kHz,  has  the  shape  of  a  cylinder  with  two  hemi¬ 
spherical  end  caps,  and  is  71  mm  long  and  41  mm  in  diameter  with  a 
mass  of  54  g.  This  sensor  is  nearly  neutrally  buoyant  in  water  so  as  to 
have  the  best  motion  coupling  with  the  water. 

Figure  1  shows  the  tank  configuration.  The  vector  sensor  was  sus¬ 
pended  near  the  center  using  a  tether  arrangement  similar  to  that 
shown  in  Shipps  and  Deng  (2003,  their  Figure  2).  This  tether  ar¬ 
rangement  had  a  resonant  frequency  of  approximately  1  Hz,  well 
outside  of  our  data  band,  and  allowed  the  sensor  to  couple  well  with 
the  local  particle  motion.  The  electrical  data  cables  were  left  slightly 
slack  to  minimize  their  interference  with  the  motion  of  the  sensor. 
The  sound  source  was  a  single  transducer  driven  using  a  synthesized 
Ricker  wavelet  with  a  median  frequency  of  7  kHz,  moved  along  sev¬ 
eral  straight  lines  while  the  sensor  remained  stationary. 

After  shooting  several  reference  lines,  two  air-filled  floats  were 
placed  in  the  water  as  highly  reflective  scattering  targets.  All  of  the 
data  presented  in  this  article  have  the  target  floats  in  the  same  hori¬ 
zontal  plane  as  the  source  and  receiver,  whereas  the  processing  and 
imaging  retained  the  3D  approach.  This  is  to  present  the  concepts  of 
vector  acoustic  data  in  the  simplest  form. 
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Figure  2  shows  common-receiver-point  gathers  of  the  accelerom¬ 
eter  and  hydrophone  data.  The  direct  arrival  from  the  source  is  the 
most  prominent  signal  and  is  hyperbolic  with  a  phase  shift  apparent 
in  the  apex  of  the  y-component  data.  The  main  pulse  of  the  signal  has 
a  delay  of  1.5  ms  from  the  time  break.  Unfortunately,  this  1.5-ms  de¬ 
lay  is  filled  with  noise  caused  by  flaws  in  the  synthesis  of  the  Ricker 
wavelet.  This  noise  degrades  the  images  shown  later.  The  sandy  bot¬ 
tom  of  the  tank  and  the  water-surface  reflections  are  visible  as  hyper¬ 
bolic  reflections,  with  their  apexes  at  the  same  offset  as  the  direct  ar¬ 
rival  and  times  of  3.9  and  4.5  ms.  The  scattered  waves  from  target 
floats  in  the  tank  show  up  as  similar  hyperbolic  reflections  with 
apexes  that  are  offset  to  the  sides  and  beginning  times  of  about 
3.5  ms. 


DIRECTIONALITY,  INTENSITY,  AND  FIDELITY 

Acoustic  intensity  I  is  the  energy  flux  caused  by  passing  sound 
waves,  quantified  by 

I  =  /?v,  (1) 

where  p  is  the  pressure  and  v  is  the  particle  velocity  (Pierce,  1981, 
equation  1-11 .3).  Ideally,  we  would  measure  the  acoustic  intensity  at 
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Figure  1 .  The  experimental  setup  was  in  a  water  tank  with  a  sand  bot¬ 
tom  designed  specifically  for  acoustic  experiments.  The  receiver 
consists  of  a  three-axis  accelerometer  and  a  hydrophone.  The  receiv¬ 
er  was  suspended  by  several  thin  elastic  strands  (not  shown)  and  was 
stationary  throughout  the  experiment.  The  acoustic  transponder 
source  was  attached  to  a  robot  that  could  move  in  three  dimensions 
and  was  computer  controlled  for  automated  data  collection.  The  tar¬ 
gets  were  air-filled  floats  and  were  placed  in  several  positions.  The 
data  in  this  article  were  taken  with  both  targets  in  the  same  horizontal 
plane  as  the  source  and  receiver.  The  receiver,  targets,  and  source 
were  2  m  below  the  water  surface;  the  three  source  lines,  shown  as 
dashed  lines,  were  offset  1  m  from  the  receiver,  the  vector  sensor  re¬ 
ceiver  was  at  (1 .527, 1.19, 0),  and  the  targets  were  at  (2.417, 1 .85, 0) 
and  (2.497, 0  64, 0)  in  (jc,y,2)  coordinates  in  meters. 


Figure  2.  Time-series  offset  plots  of  the  accelerometer  components 
and  the  hydrophone  for  a  data  set  acquired  while  moving  the  source 
along  the  y-axis  from  (0.52,  — 1.1 1,0)  to  (0.52, 2.29,0)  in  (jc,y,z)  co¬ 
ordinates  in  meters.  The  direct  arrival  from  the  source  appears  as  the 
prominent  hyperbolic  signal  in  all  channels  with  the  y-axis  node  and 
phase  shift  at  the  hyperbola’s  apex.  The  main  pulse  of  the  signal  has  a 
delay  of  1 .5  ms  from  the  time  break.  The  sandy  bottom  of  the  tank 
and  the  water-surface  reflections  also  appear  as  hyperbolic  reflec¬ 
tions,  particularly  in  the  2-axis  panel,  with  their  apexes  at  the  same 
offset  as  the  direct  arrival  and  traveltimes  of  3.9  and  4.5  ms.  The 
scattered  waves  from  the  target  floats  appear  as  hyperbolic  reflec¬ 
tions  with  apexes  at  y  =  0.64  and  1.85  m  and  traveltimes  of  about 
3.5  ms.  The  scattered  waves  from  the  targets  are  about  1  %  of  the  am¬ 
plitude  of  the  direct  wave. 


3D  vectoracoustic  imaging 


Q3 


every  sensor  location,  but  there  is  no  such  ideal  sensor.  Motion  sen¬ 
sors,  such  as  the  accelerometers  used  in  this  study,  generally  have 
vector  fidelity  problems,  and  hydrophones  usually  measure  the  ve¬ 
locities  of  their  faceplates  rather  than  the  true  pressure. 

When  imaging,  we  do  not  need  to  calculate  the  true  intensity,  yet 
we  need  some  value  that  resembles  it  in  terms  of  direction,  phase, 
and  amplitude.  For  a  single-axis  pressure  gradient  sensor,  the  inten¬ 
sity  component  is  defined  as 

In  =  P1P2  sin(0!  -  (t>2)/2p0cod ,  (2) 

where  Px  and  P2  are  the  pressures  at  the  transducer  pairs,  <f>]  and  <f>2 
are  the  signal  phases  at  the  transducer  pairs,  p0  is  the  mean  density,  tv 
is  the  angular  frequency,  and  d  is  the  transducer  spacing  (Fahy,  1 989, 
equation  5.7). 

For  accelerometer  sensors,  there  is  no  precise  definition  of  inten¬ 
sity.  Directionality  (the  accuracy  of  the  direction  estimate)  depends 
on  good  matching  of  the  phase  and  amplitudes  of  each  axis  of  the  ac¬ 
celerometers  or  of  each  hydrophone  element  of  a  pressure  gradient 
sensor.  The  directionality  is  degraded  by  mismatched  phases,  ampli¬ 
tude  sensitivity,  environmental  noise,  and  vector  infidelity.  Direc¬ 
tionality  also  can  be  lost  when  two  waves  from  different  directions 
arrive  simultaneously. 

Figure  3  shows  a  hodograph  of  the  direct  arrivals  from  the  x  and  y 
accelerometer  components  with  the  source  at  four  different  loca¬ 
tions.  This  figure  illustrates  some  of  the  inaccuracies  of  the  sensor 
and  of  the  experimental  setup.  Two  of  the  signals  are  aligned  with  the 
x -  and  y-axes  and  two  are  at  45°  angles  from  the  axes.  The  y-axis  ac¬ 
celerometer  of  this  sensor  is  about  twice  as  sensitive  as  the  jc-axis  ac¬ 
celerometer,  and  the  data  shown  arc  corrected  for  this.  The  sensor 
motion  from  the  source  in  the  y-axis  direction  (green)  was  aligned 
well  with  the  y-axis,  but  the  sensor  motion  from  the  source  in  the  x 
direction  (red)  deviated  substantially  from  the  axis.  The  x-direction 


Figure  3,  Hodograph  of  the  x  and  y  accelerometer  components  of  di¬ 
rect  waves  from  the  source.  Two  of  the  signals  are  from  the  source 
aligned  with  the  x-  and  y-axis,  and  two  are  from  the  source  at  45°  an¬ 
gles  from  the  axes.  The  y-axis  accelerometer  was  twice  as  sensitive 
as  the  jc-axis  accelerometer,  and  the  data  have  been  corrected  for  this 
difference. 


motion  deviations  may  have  been  caused  by  an  internal  defect  in  the 
sensor  or  by  some  motion  restrictions  from  the  electrical  cables.  In 
spite  of  the  vector  fidelity  problems  with  our  sensor,  we  are  able  to 
use  these  data  for  this  demonstration  of  the  vector  acoustic  imaging 
concept. 


DETECTION  AND  LOCATION 


The  location  of  a  controlled  source  can  be  determined  with  a  sin¬ 
gle  vector  sensor  by  measuring  the  direction  and  the  traveltime  of  the 
direct  wave.  The  source  position  vector  s  is 


s  =  g  +  d  tc,  (3) 

where  g  is  the  receiver  position,  d  is  the  unit  vector  aligned  with  the 
particle  motion,  t  is  the  traveltime,  and  c  is  the  sound  speed. 

Figure  4  shows  how  a  scattering  target  can  be  located  using 
known  source  and  receiver  locations  by  measuring  the  directional 
response  and  the  traveltime  of  the  scattered  waves.  Multiple  scatter¬ 
ing  targets  can  be  located  with  a  single  vector  sensor,  provided  that 
the  scattered  signals  are  separated  in  time  at  the  sensor.  In  the  follow¬ 
ing  calculations,  the  source  and  receiver  positions  s  and  g  are  known, 
and  the  traveltime  t  and  the  vector  response  of  the  sensor  d  are  mea¬ 
sured.  The  target  position  vector  e  is  calculated  by 


and 


e  =  g  +  dr, 

b~  -  c  r 


2  (b  cos  0  -  ct)  ’ 


(4) 

(5) 


0  =  cos 


d  •  (g  —  s) 
b\d\ 


(6) 


where  r  is  the  range  to  the  target,  b  is  the  source  to  receiver  range,  d  is 
the  data  vector,  d  is  the  data  vector  of  unit  length,  6  is  the  angle  be¬ 
tween  the  receiver-source  and  the  receiver-target  vectors,  and  II  de¬ 
notes  the  norm,  or  length,  of  the  vector.  The  large  dot  denotes  the 
vector  dot  product. 

Figure  5  shows  the  sensor  motions  caused  by  the  scattered  waves 
from  the  targets.  This  hodograph  illustrates  the  directional  informa¬ 
tion  in  vector  acoustic  data.  Interference  from  the  direct  signals  is 
present,  even  though  these  signals  were  selected  for  minimum  inter¬ 
ference.  The  interference  from  the  direct  signal  will  cause  loss  of  res¬ 
olution  in  the  target  imaging  and  also  may  be  confused  with  poor 
vector  fidelity. 


Figure  4.  Geometry  of  the  target  location  The  source  and  sensor  lo¬ 
cations  s  and  g  are  known,  along  with  their  separation  b.  The  record¬ 
ed  data  d  are  the  vector  responses  of  the  sensor. 
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When  multiple  scattered  signals  overlap,  such  as  in  the  data  in 
Figure  2,  then  the  targets  and  surfaces  cannot  be  located  with  a  single 
source  and  receiver  location.  Locating  each  target  within  a  field  of 
scattering  targets  requires  multiple  sources  or  sensors.  For  most  cas¬ 
es,  a  single  line  of  either  source  or  vector  receiver  locations  is  suffi¬ 
cient.  A  single  sensor  or  receiver  will  not  be  able  to  correctly  locate 
two  scattering  targets  that  are  positioned  symmetrically  with  respect 
to  the  axis  of  the  line.  This  two-target  ambiguity  can  be  resolved  in 
several  ways.  If  either  the  source  or  the  vector-sensor  is  held  fixed  as 
the  other  is  moved  along  a  line,  the  fixed  device  can  be  off-axis  of  the 


Figure  5.  Hodograph  of  the  x  and  y  accelerometer  components  of  the 
scattered  signals  from  the  target  floats.  The  traces  were  selected  at  a 
location  where  the  reflections  had  the  best  separation  from  other  sig¬ 
nals.  The  scattered  arrivals  were  time  windowed  and  corrected  for 
the  different  sensitivity  of  the  x-  and  y-axis  accelerometers.  There  is 
some  remaining  energy  from  the  direct  signals  that  interferes  with 
these  scattered  signals  and  shows  up  as  the  weaker,  sideways  wig¬ 
gling  of  the  signals.  The  source  positions  were  (  —  1.0,  -0./5)  and 
( - 1 .0,0.9),  and  the  sensor  was  at  (0.0, 0.0). 


Figure  6.  Time  series  plots  of  the  three  accelerometer  signals  and  the 
hydrophone.  The  hydrophone  response  must  be  corrected  to  match 
the  response  of  the  accelerometers  before  using  the  hydrophone 
trace  to  rectify  the  accelerometer  traces. 


line.  This  reduces  the  location  ambiguity  to  objects  at  symmetric  po¬ 
sitions  on  a  single  plane. 

All  target  position  ambiguities  can  be  eliminated  by  moving  either 
the  source  or  the  receiver  (or  both)  along  a  serpentine  path.  With  a 
serpentine  path,  there  never  will  be  two  target  positions  that  always 
have  identical  two-way  traveltimes.  This  presumes  that  the  positions 
of  the  source  and  the  receiver  and  the  shot  times  are  always  known.  A 
vector  sensor  can,  in  fact,  be  used  to  position  itself  relative  to  a  sound 
source  if  it  includes  orientation  devices,  such  as  a  compass  and  tilt- 
meters.  This  eliminates  the  need  for  additional  positioning  systems 
that  are  necessary  for  3D  surveys. 

Rectifying  the  signal 

The  direction  vector  d  in  equation  3  is  proportional  to  the  norm  of 
the  time-averaged  acoustic  intensity  I  (D’ Spain  et  al.,  1 99 1 ;  Crocker 
and  Arenas,  2003)  or  to  the  direction  of  energy  flux  vector  (Chap¬ 
man,  2004)  and  is  analogous  to  the  Poynting  vector  used  in  electro¬ 
magnetism.  Although  the  acoustic  intensity  would  be  ideal  for  imag¬ 
ing,  the  actual  sensor  measurements  are  the  instantaneous  particle 
accelerations  and  oscillate  between  the  positive  and  negative  inten¬ 
sity  vectors. 

The  instantaneous  accelerations  can  be  easily  rectified  to  give  ap¬ 
proximate  intensity  vectors  without  excessive  calculations  by  using 
the  hydrophone  data  in  conjunction  with  the  accelerometer  data.  The 
sensor  we  used  has  an  integral  hydrophone,  but  its  response  has  a  dif¬ 
ferent  phase  and  different  frequency  than  the  accelerometers  (Figure 
6).  The  hydrophone  signal  must  first  be  filtered  before  it  can  be  used 
to  rectify  the  accelerometer  signals,  as  quantified  by 


a(t)  *-*  A{J) 

(7) 

(8) 

where  a(t)  and  p(t)  are  the  time  series  signals  of  the  time-windowed 
direct  arrivals  from  the  accelerometer  and  the  pressure  hydrophone 
and  where  A{f)  and  P(f)  are  their  respective  Fourier  transforms  into 
the  frequency  domain. 

The  original  hydrophone  spectrum  PD(f)  is  whitened  by  dividing 
it  by  the  direct  signal  spectrum  P(f),  then  multiplied  by  the  acceler- 


Time  (ms) 


Figure  7.  Time  series  plot  of  the  y-axis  accelerometer,  the  pressure 
hydrophone,  and  the  spectrum-corrected  (filtered)  hydrophone  sig¬ 
nal  for  the  direct  arrival  only.  The  filtered  hydrophone  data  are  plot¬ 
ted  slightly  above  the  accelerometer  data.  There  is  an  almost  exact 
match  between  the  accelerometer  waveform  and  the  filtered  pres¬ 
sure  waveform. 
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ometer  direct  signal  spectrum  A(f)  to  give  a  filtered  hydrophone 
spectrum  PA{f)  that  can  be  transformed  to  the  amplitude  and  phase- 
matched  pressure  time  series pA(t)  shown  in  Figure  7.  This  sequence 
is  represented  by 


PM)  -  W  • 

(9) 

PA(f)  =  Pvtf-AU), 

(10) 

PA(t)*~*  PaU)- 

(11) 

The  matched  pressure  time  series  is  used  to  create  the  rectifying 
sign  series  sgn(pA(/))  that  can  change  the  accelerometer  data  aD(i) 
into  rectified  accelerometer  data  aR(l ),  as  in 

aR(t)  =  sgn[/?A(f)]  aD(t).  (12) 

The  three  rectified  accelerometer  components  can  then  be  com¬ 
bined  in  the  acoustic  intensity  vector  I  (equation  1 ).  If  no  pressure 
data  are  available,  the  accelerometer  data  can  be  partially  rectified,  if 
the  source  wavelet  is  known,  by  using  the  crosscorrelation  of  the 
source  wavelet  a{i)  and  the  accelerometer  data  aD(t),  given  by 

aR(t)  =  sgn[aD(r)  *  a(t)]  •  aD(t) .  (13) 

As  illustrated  in  the  lower  plot  in  Figure  8,  equation  13  will  not 
work  for  phase-shifted  reflections,  such  as  from  the  water  surface  or 
air  bubbles. 

IMAGING 

Each  source  and  receiver  location  s,  g  produces  a  3-C  trace, 
d(s,g,f).  Equations  3  and  4  map  each  sample  from  d(s,g,/)  into  a  tar¬ 
get  position  vector  e(s,g,/).  A  volume  image  can  be  built  by  sum¬ 
ming  the  absolute  amplitudes  |d(s,g,/)|  for  all  of  the  recorded  traces 
into  their  respective  volume  bins. 

Figure  9  shows  horizontal  slices  through  three  volume  images 
made  from  three  different  source  lines  at  the  sensor  depth,  showing 
the  two  float  targets  as  well  as  more  distant  walls.  The  direct  arrivals 
have  been  muted  before  imaging.  Most  of  the  targets  and  walls  have 
ghosts  from  the  unrectified  negative  parts  of  the  signals  The  noise 
shows  up  as  uncertainty  in  the  direction  to  the  targets  but  not  in  the 
ranges.  The  data  were  corrected  for  different  sensitivities  of  the  sen¬ 
sor  components  before  imaging. 

The  ghosts  and  some  of  the  noise  present  in  Figure  9  can  be  re¬ 
duced  by  rectifying  the  data  using  the  hydrophone  data  as  described 
by  equations  9-1 2  and  illustrated  in  Figures  7  and  8.  Figure  10  illus¬ 
trates  the  imaging  in  the  same  horizontal  slices  as  in  Figure  9.  The 
ghosts  and  the  noise  are  much  weaker,  and  the  targets  and  walls  are 
positioned  correctly. 

DISCUSSION 

This  article  demonstrates  a  method  for  acoustically  imaging  a 
very  simple  environment.  A  real  marine  environment  is  far  more 
complex,  with  many  surfaces  and  objects  that  pose  challenges  to  the 
processing  of  vector  acoustic  data.  There  will  need  to  be  a  great  deal 
of  theoretical  development  in  vector  data  processing  and  imaging 
before  this  method  can  be  applied  to  realistic  environments.  The  pri¬ 


mary  challenges  are  the  multiple  overlapping  reflections  and  the  un¬ 
known  sound  velocity  structure.  Vector  wavefield  decomposition 
methods  can  separate  and  isolate  the  individual  reflections  (e.g., 
Muijs  et  al.,  2004;  Paulus  and  Mars,  2006)  that  can  then  be  imaged. 
Alternatively,  existing  migration  concepts  can  be  extended  to  vector 
data  (e.g.,  Liith  et  al.,  2005).  It  also  may  be  possible  to  determine  ve¬ 
locities  using  backscattered  energy  (Wood  et  al .,  1 997)  that  could  al¬ 
low  for  surveys  with  fewer  receiver  locations  and  smaller  spatial  ap¬ 
ertures. 

The  direct  imaging  method  shown  in  this  article  does  not  use  any 
phase  information  between  source  positions,  and  the  images  do  not 
have  some  of  the  common  artifacts  because  of  frequency  and  spatial 
limitations  of  the  data.  The  primary  flaws  in  these  vector-derived  im¬ 
ages  are  the  directional  noise  caused  by  overlapping  signals  and  the 
limited  illumination  of  the  flat  surfaces.  Incorporating  vector  data 
into  existing  or  future  imaging  methods  may  reduce  artifacts  signifi¬ 
cantly,  particularly  when  the  spatial  coverage  is  limited. 

The  illumination  of  surfaces  and  extended  structures  for  vector- 
based  surveys  is  an  unexplored  issue.  Low-frequency  seismic  sur¬ 
veys  over  sedimentary  basins  produce  primarily  specular  reflections 
and  may  require  more  extensive  source  and  receiver  coverage.  Some 


Figure  8.  Time  series  plots  of  the  z- axis  accelerometer,  the  pressure 
hydrophone,  and  the  rectified  z-axis  signal,  (a)  Plot  uses  the  matched 
hydrophone  data  as  the  rectifying  series,  (b)  Plot  uses  the  wavelet 
correlated  with  the  ^-axis  data  as  the  rectifying  series.  The  (a)  pres¬ 
sure-rectified  data,  although  not  perfectly  phase  matched,  show 
most  of  the  tank  bottom  and  water  surface  reflections  coming  from 
opposite  directions.  The  (b)  wavelet-rectified  data  erroneously  place 
the  water  surface  reflection  because  of  the  77-radians  phase  shift  of 
the  free  surface  reflection. 
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higher-frequency  surveys  for  near-seafloor  structures  and  objects 
use  primarily  diffuse  (backscatter)  reflections  and  point  diffractions 
and  may  require  only  a  few  widely  spaced  survey  lines. 

Figure  1 1  compares  the  vector-imaged  volume  with  a  2D-migrat- 
ed  image  generated  from  the  hydrophone  data  shown  in  Figure  2. 


x-line  right 


The  vector  image  correctly  positions  the  scattering  points  and  re¬ 
flecting  surfaces  in  3D  space,  whereas  2D  migration  cannot  correct¬ 
ly  position  scattering  points  that  are  outside  the  horizontal  plane  of 
the  source  line  and  receiver.  This  is  not  a  fair  test  for  2D  migration 
because  the  survey  configuration  and  the  environment  in  this  experi¬ 
ment  are  not  restricted  to  the  vertical  plane  that  2D  migration  was  de¬ 
signed  for.  Also,  2D  migration  is  restricted  to  one  direction  from  the 
source  or  receiver  line,  usually  below,  whereas  the  tank  environment 
used  here  had  features  in  all  directions. 
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Figure  9.  Horizontal  slices  of  3D  volume  images  of  the  water  tank 
using  unrectified  data  from  three  different  source  lines.  The  source 
(red  lines),  sensor  (red  dots  in  center),  and  targets  (black  circles)  are 
all  at  the  same  depth  in  these  slices  in  the  x-y  plane.  All  of  the  reflect¬ 
ing  and  scattering  structures  have  strong  ghosts  because  the  data  are 
unrectified.  The  images  were  calculated  by  mapping  each  time  sam¬ 
ple  to  a  gridded  position  using  equation  4  and  adding  the  absolute 
value  of  the  amplitude  of  that  time  sample  to  the  corresponding  grid. 
The  direct  arrivals  are  not  included  in  these  images. 


Figure  10.  Rectified  (pressure-derived)  horizontal  slices  of  the  vol¬ 
ume  image  of  the  tank  environment.  The  source  (red  lines),  sensor 
(red  dots  in  center),  and  targets  (black  circles)  are  all  at  the  same 
depth  in  these  slices  in  the  x-y  plane.  The  images  were  formed  in  the 
same  manner  as  in  Figure  9.  Although  the  rectification  is  imperfect, 
most  of  the  ghosts  have  been  eliminated,  and  the  noise  has  been  re¬ 
duced  substantially. 
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waves  at  relatively  few  locations,  it  is  possible  to  image  the  acoustic 
features  of  a  3D  volume. 

The  technique  presented  here  is  limited  to  simple  environments 
with  a  small  number  of  surfaces  and  targets.  With  the  development 
of  methods  for  separating  the  overlapping  vector  signals  and  invert¬ 
ing  the  decomposed  wavefield,  it  may  be  possible  to  image  geologi¬ 
cally  realistic  3D  environments  with  significantly  reduced  source 
and  receiver  coverage.  A  3D  seismic  survey  may  be  possible  with 
vector  sensors  using  logistics  similar  to  2D  surveys  that  use  conven¬ 
tional  hydrophones. 
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